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Observation of a New X b State in Radiative Transitions to T(1S) and T(2S) at ATLAS 
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(Dated: December 21, 2011) 

The Xb(nP) quarkonium states are produced in proton-proton collisions at the Large Hadron Collider 
(LHC) at y/s — 7 TeV and recorded by the ATLAS detector. Using a data sample corresponding to 
an integrated luminosity of 4.4 fb _1 , these states are reconstructed through their radiative decays 
to T(1S,2S) with T — > fi~ . In addition to the mass peaks corresponding to the decay modes 
X 6 (1P,2P) -> T(1S) 7, a new structure centered at a mass of 10. 530±0.005 (stat.)±0.009 (syst.) GeV 
is also observed, in both the T(1S) 7 and T(2S) 7 decay modes. This structure is interpreted as the 
Xb(3P) system. 

PACS numbers: 12.38.-t, 13.20.Gd, 14.40.Pq, 14.65.Fy 



Measurements of the properties of heavy quark- 
antiquark bound states, or quarkonia, provide a unique 
insight into the nature of Quantum Chromodynamics 
close to the strong decay threshold. For the bb sys- 
tem, the quarkonium states with parallel quark spins 
(s = 1) include the S-wave T and the P-wave X b states, 
where the latter each comprise a closely spaced triplet 
of J = 0, 1,2 spin states: X b o, "Xbi and Xbi- The Xfe(LP) 
and Xb(2P), with spin- weighted mass barycenters of 9.90 
and 10.26 GeV, respectively, can be readily produced in 
the radiative decays of T(2S) and T(3S) and have been 
studied experimentally 

In this letter, Xb quarkonium states are reconstructed 
with the ATLAS detector through the radiative decay 
modes X b (riP) -> T(1S)7 and X b (nP) -> T(2S)7, in 
which T(1S, 2S) — > and the photon is recon- 

structed either through conversion to e + e~ or by direct 
calorimetric measurement. Previous experiments have 
measured the X{,(1P) and Xf,(2P) through these decay 
modes |2| . The X b (3P) state has not previously been ob- 
served. It is predicted to have an average mass of approx- 
imately 10.52 GeV, with hyperfine mass splitting between 
the triplet states of 10-20 MeV @,|||. 

The ATLAS detector 5] is a general-purpose parti- 
cle physics detector with a forward-backward symmetric 
cylindrical geometry and near 47r coverage in solid an- 
gle. The inner tracking detector (ID) consists of a silicon 
pixel detector, a silicon microstrip detector, and a transi- 
tion radiation tracker (TRT) . The ID is surrounded by a 
thin superconducting solenoid providing a 2 T magnetic 
field and by high-granularity liquid-argon sampling elec- 
tromagnetic (EM) calorimeters. An iron-scintillator tile 
calorimeter provides hadronic coverage in the central ra- 
pidity range. The endcap and forward regions are instru- 
mented with liquid-argon calorimeters for both electro- 
magnetic and hadronic measurements. The muon spec- 
trometer (MS) surrounds the calorimeters and consists 
of a system of precision tracking chambers and detectors 
for triggering, inside a toroidal magnetic field. 

The data sample used for this measurement was 



recorded by the ATLAS experiment during the 2011 LHC 
proton-proton collision run at a center-of-mass energy 
y/s = 7 TeV. The integrated luminosity of the data sam- 
ple, which includes only data-taking periods where all rel- 
evant detector sub-systems were operational, is 4.4 fb -1 . 
A set of muon triggers designed to select events con- 
taining muon pairs or single high transverse momentum 
muons was used to collect the data sample. 

In this analysis each muon candidate must satisfy stan- 
dard muon quality requirements [f|. It must have a 
track, reconstructed in the MS, combined with a track 
reconstructed in the ID with transverse momentum px > 
4 GeV and pseudorapidity |ry| < 2.3. The di-muon selec- 
tion requires a pair of oppositely charged muons, which 
arc fitted to a common vertex. A very loose vertex qual- 
ity requirement (x 2 per degree of freedom [d.o.f.] < 20) 
is used and no mass or momentum constraints are ap- 
plied to the fit. The di-muon candidate is also required 
to have pt > 12 GeV and \y\ < 2.0. The invariant 
mass distribution, m^, of di-muon candidates is shown 
in Fig. [T] Those candidates with masses in the ranges 
9.25 < m w < 9.65 GeV and 9.80 < m MM < 10.10 GeV 
are selected as T(1S) — > H + p~ an d T(2S) — > [i + p~ can- 
didates respectively. The asymmetric mass window (evi- 
dent from Fig.[IJ for T(2S) candidates is chosen in order 
to reduce contamination from the T(3S) peak and con- 
tinuum background contributions. 

The reconstruction of photons in ATLAS is described 
in Ref. Q . Further details related to this particular anal- 
ysis are described below. 

Converted photons are reconstructed from two oppo- 
sitely charged ID tracks intersecting at a conversion ver- 
tex, with the opening angle between the two tracks at this 
vertex constrained to be zero. For tracks with signals in 
the TRT, the transition radiation should be consistent 
with an electron hypothesis. In order to be reliably re- 
constructed, each conversion electron track must have a 
minimum transverse momentum of 500 MeV. It is also 
required to have at least four silicon detector hits and not 
to be associated to either of the two muon candidates. To 
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FIG. 1: The invariant mass of selected di-muon candidates. 
The shaded regions A and B show the selections for T(1S) 
and T(2S) candidates respectively. 



reduce background contamination, the conversion candi- 
date vertex is required to be at least 40 mm from the 
beam axis and have a vertex \ 2 probability of greater 
than 0.01. The converted photon impact parameter with 
respect to the di-muon vertex is required to be less than 
2 mm. 

Electromagnetic calorimeter energy deposits not 
matched to any track are classified as unconverted pho- 
tons. This analysis uses the "loose" photon selection de- 
scribed in Ref. [7J, with a minimum photon transverse 
energy of 2.5 GeV. The loose photon selection includes a 
limit on the fraction of the energy deposit in the hadronic 
calorimeter as well as a requirement that the transverse 
width of the shower be consistent with the narrow shape 
expected for an EM shower. 

To check that an unconverted photon originates from 
the same vertex as the T, and to improve the mass reso- 
lution of the reconstructed Xb, the polar angle of the pho- 
ton is corrected using the procedure described in Ref. [§[ . 
The corrected polar angle is determined using the mea- 
surement of the photon direction from the longitudinal 
segmentation of the calorimeter and the constraint from 
the di-muon vertex position. Photons incompatible with 
having originated from the di-muon vertex are rejected 
by means of a loose cut on the fit result (x 2 per d.o.f. 
< 200). 

The converted (unconverted) photon candidates are re- 
quired to be within \i]\ < 2.30(2.37). Unconverted pho- 
tons must also be outside the transition region between 
the barrel and the endcap calorimeters, 1.37 < < 1.52. 

The Xb candidates are formed by associating a recon- 
structed T — > candidate with a reconstructed pho- 
ton. The invariant mass difference Am = m(fi + fi~ 7) — 
m{jjL + n~) is calculated to minimize the effect of T — > 



mass resolution. In order to compare the Am 
distributions of both X b (nP) ->■ T(1S)7 and X b (nP) 
T(2S)7 decays, the variable m^ = Am + mx(fes) is de- 
fined, where mj-(ks) are the world average masses [t| of 
the Y(fcS) states. Requirements of pr(M + M~) > 20 GeV 
and pt(m + aO > 12 GeV are applied to T candidates 
with unconverted and converted photon candidates re- 
spectively. These thresholds are chosen in order to opti- 
mize signal significance in the Xf,(lP,2P) peaks. 

Figure [2f a) shows the mi distribution for unconverted 
photons and Fig. [2fb) the mi and rhi distributions for 
converted photons. In addition to the expected peaks for 
Xb(lP,2P) —> T(1S,2S)7, structures are observed at an 
invariant mass of approximately 10.5 GeV. These addi- 
tional structures are interpreted as the radiative decays 
of the previously unobserved X&(3P) states, X&(3P) — > 
T(1S) 7 and X h (3P) ->• T(2S) 7 . 

Separate fits are performed to the m^ distributions of 
the selected candidates reconstructed from con- 

verted and unconverted photons to extract mass informa- 
tion from the observed Xb(3P) signals. The higher thresh- 
old for unconverted photons (2.5 GeV, versus 1 GeV for 
converted photons) prevents the reconstruction of the 
soft photons from X(,(2P,3P) decays into T(2S). 

An unbinned extended maximum likelihood fit is per- 
formed to the mi = Am + mx(is) distribution of the se- 
lected unconverted /i + /i~7 candidates. The three peaks 
in the distribution are each modeled by a Gaussian prob- 
ability density function (pdf) with independent normal- 
ization parameter, N n , mean value, m„, and width pa- 
rameter, er„. The background distribution is parameter- 
ized by the pdf Ve-exp (A • Am + B ■ Am -2 ) where Nb, 
A, and B are all free parameters. The three mean values 
mn=i,2,3 determined by the fit are shown in TablelU The 
mean value 7713 is an estimate of the mass barycenter of 
the observed Xb(3P) signal. 

Likewise, the mi = Am + my(is) and m.2 = Am + 
TO T(2S) distributions for the sample of /i + // - 7 candidates 
reconstructed from converted photons are fitted using an 
unbinned extended maximum likelihood method. A si- 
multaneous fit is performed on the mi and m.2 distri- 
butions for the X b {nP) ->■ T(1S)7 (for n = 1,2,3) and 
X b {nP) -> T(2S)7 (for n = 2,3 only) signals, with the 
distributions modeled by three signal components (two 
of which are shared between the T(1S) and T(2S) distri- 
butions) and two background distributions. 

In the Am distribution for the converted photon can- 
didates the typical mass resolution is found to be in the 
range 16 — 20 MeV, of similar magnitude to the hypcrfine 
splittings, motivating the need for multiple signal com- 
ponents for each of the Xb{nP) peaks. For n = 1,2, the 
radiative branching fractions of the J = states are sup- 
pressed with respect to the J = 1, 2 states Q and there- 
fore a J = component is not included in the fit. Similar 
behavior is assumed for the n = 3 case. Each of the three 
peaks (n — 1,2,3) is therefore parameterized by a dou- 
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TABLE I: The fitted mass of the Xb(nP) signals for both converted and unconverted photons. The systematic uncertainty 
on the mass of candidates reconstructed with unconverted photons is determined in the same way for all three states. Also 
included are theoretical predictions [3, 13 f° r the spin-averaged masses of the Xb states. 



State 


Model predictions 
[3, 4] [MeV] 


Unconverted Photons 


Fitted masses 


[MeV] 

Converted Photons 




Xt(lP) 
X b (2P) 
X b {3P) 


9900 
10260 
10525 


9910 ± 6 (stat.) ± 11 (syst.) 
10246 ± 5 (stat.) ± 18 (syst.) 
10541 ± 11 (stat.) ± 30 (syst.) 


Fixed to 
Fixed to 
10530 ± 


Xbi = 9892.78 & X b2 = 9912.21 
Xbi = 10255.46 & X b2 = 10268.65 
5 (stat.) ± 9 (syst.) 


[9] 
[9] 



blet of Crystal Ball (CB) [l(| functions (corresponding 
to J = 1 , 2 states) with resolution a and radiative tail 
parameters common to all peaks. For n ~ 1 and n = 2, 
the peak mass values and hyperfine splittings are fixed 
to the world averages Q for the respective Xb states (see 
Table Q} . For n = 3, the hyperfine mass splitting is fixed 
to the theoretically predicted value of 12 MeV [4j, while 
the average mass is left as a free parameter. The un- 
known relative normalization of the J = 1 and J = 2 CB 
peaks is taken to be equal and treated as a systematic 
uncertainty (for all doublets) for the baseline fit. 

In order to take into account energy loss from the 
photon conversion electrons due to bremsstrahlung and 
other processes, the measured values of Am in the fh\ 
and 77i2 distributions are scaled by a common param- 
eter A = 0.961 ± 0.003, which determines the energy 
scale and is derived from the fit to the Xb(lP, 2P) signals. 
The background components of the Am distributions for 
the T(1S)7 and T(2S)7 final states are each modeled 
by the pdf N% ■ (Am - q° k ) Ak ■ exp [B k ■ (Am - q° k )] for 
Am > , and zero otherwise, where Ng , g° , , and Bk 
(k = 1,2) are all free parameters. The mean value W13 
determined by the fit is shown in Table HI 

In the fit using unconverted photons, the signal is refit- 
ted using an alternative (two Gaussians) model for each 
of the three Xb states, resulting in a negligible change 
in the peak positions. Alternative fits to the back- 
ground are also used, either including constraints on the 
Am distribution using di-muon pairs from the low-mass 
(8.0 GcV < m w < 8.8 GeV) sideband or different back- 
ground pdfs. The systematic uncertainty on the Xb(3P) 
mass barycenter from modeling of the background distri- 
bution is determined to be ±21 MeV. The systematic un- 
certainty associated with the unconverted photon energy 
scale is estimated to be ±2% on the Am position, corre- 
sponding to a systematic uncertainty on of ±22 MeV. 
The uncertainties due to background modeling and pho- 
ton energy scale comprise the dominant sources of sys- 
tematic uncertainty. 

For the fit using converted photons, alternative signal 
and background models are compared, as well as releas- 
ing various constraints in the fit model. The unknown 
relative normalizations of the J — 1 and J — 2 CB peaks 
are varied both coherently and incoherently between the 
IP, 2P and 3P doublets by ±0.25, resulting in a max- 



imum variation in 7713 of ±5 MeV. Smaller variations 
are obtained if the common value of the relative nor- 
malization is allowed to be determined freely by the fit 
to the three doublets. Background modeling variations, 
decoupled fits to the rh\ and m 2 distributions, and indi- 
vidually released constraints on the mass position of the 
n = 1,2 doublets each result in deviations of the order 
of ±5 MeV or smaller. Furthermore, if the constraints 
on the masses of the n = 1,2 peaks are released, the 
values obtained from the fit are consistent with expecta- 
tions 0, within statistical errors and uncertainty in the 
relative contributions from J = 1 and J = 2 states. The 
effect of symmetrizing the T(2S) mass window is stud- 
ied and found to have negligible effect on the fitted Xb 
masses while increasing background contamination. The 
resulting shifts in 77T3 for these independent variations 
are added in quadrature to provide an estimate of the 
systematic uncertainty. 

The Xb(3P) signal significance is assessed from 
log(L max /Lo), where L miiX and Lq are the likelihood val- 
ues from the nominal fit and from a fit with no Xb(3P) 
signal included, respectively. The fit is repeated with 
each of the systematic variations in the model, as dis- 
cussed above, and the likelihood ratio re-evaluated. The 
significance of the Xf,(3P) signal is found to be in excess 
of six standard deviations in each of the unconverted and 
converted photon selections independently. 

The mass barycenter for the Xb(3P) signal, determined 
from the fit using unconverted photon candidates is: 

m 3 = 10.541 ± 0.011 (stat.) ± 0.030 (syst.) GeV. 

The mass barycenter for the Xb(3P) signal, determined 
from the fit using converted photon candidates is: 

m 3 = 10.530 ± 0.005 (stat.) ± 0.009 (syst.) GeV. 

The measured mass barycenters of the Xt(lP), Xb(2P) 
and Xf,(3P) systems are summarized in Table HI The re- 
sults of the converted and unconverted photon analyses 
for the X&(3P) are found to be compatible. Given the 
substantially smaller systematic uncertainties in the con- 
version measurement, the final mass determination for 
S3 is quoted solely on the basis of this analysis. 

In conclusion, the production of the heavy quarkonium 
states Xb(nP) in proton-proton collisions at = 7 TeV 
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FIG. 2: (a) The mass distribution of Xb T(1S)7 candi- 
dates for unconverted photons reconstructed from energy de- 
posits in the electromagnetic calorimeter (xlt/d.o.f. = 0.85). 
(b) The mass distributions of Xb — > T(kS)y (k = 1,2) candi- 
dates formed using photons which have converted and been 
reconstructed in the ID (xlt/d.o.f. = 1.3). Data are shown 
before the correction for the energy loss from the photon con- 
version electrons due to bremsstrahlung and other processes. 
The data for decays of Xb — > T(1S)7 and Xb — > T(2S)7 are 
plotted using circles and triangles respectively. Solid lines rep- 
resent the total fit result for each mass window. The dashed 
lines represent the background components only. 



is observed through reconstruction of the radiative decay 
modes of Xb(nP) — > T(1S,2S) 7. Mass peaks correspond- 
ing to Xh(lP,2P) decays are observed, together with ad- 
ditional structures at higher mass, which are consistent 
with theoretical predictions for Xb(3P) — > T(1S)7 and 
Xf,(3P) — > T(2S)7. These observations are interpreted 
as the X&(3P) multiplet, the mass barycenter of which is 
measured to be 10.530±0.005 (stat.) ±0.009 (syst.) GeV. 
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